The first equation expresses the mass balance in the annular passage. The second and third equations represent the momentum balance in the meridional and tangential directions, respectively. The last equation represents thermal energy balance, incorporating viscous dissipation and wall heat transfer. The viscous shear stress is modeled using a conventional friction factor approach:
τ=21cfρv2
where cf is the skin friction coefficient, obtained from empirical correlations as a function of the Reynolds number. Thermodynamic properties are evaluated from the equation of state, using pressure–density as input pair.
"R. Agromayor, B. Müller, and L. O. Nord, “One-dimensional annular diffuser model for preliminary turbomachinery design,” International Journal of Turbomachinery, Propulsion and Power, vol. 4, no. 3, p. 31, 2019, doi: 10.3390/ijtpp4030031.
New model
An equivalent formulation of the model can be derived from the total energy equation, leading to the following system of governing equations:
To close the system, the density is expressed as a function of pressure and enthalpy:
dmdρ=(∂p∂ρ)hdmdp+vm(∂h∂ρ)pdmdh
where the partial derivatives of density can be expressed in terms of the speed of sound and the Gruneisen parameter:
(∂h∂ρ)p=−c2ρG(∂p∂ρ)h=c21+G
Substituting the above thermodynamic relations into the continuity equation gives the following quasi-linear system of ordinary differential equations:
This result shows that the system becomes singular when the meridional Mach number equals unity. At this condition, the flow reaches a sonic state along the mean streamline, corresponding to choking at the section of minimum area.
Total pressure and entropy generation
The mechanical energy equation with this term is given by:
This equation makes a lot of sense. The decrease in stagnation pressure is caused by the wall friction, plus any other sink terms that we put in the momentum equations. In this case we see clearly that I has the physical meaning of stagnation pressure drop.
Their effect on the entropy generation expression is given by
σT=b2(1−TwT)qw+b2τv+Ivm
In general, which even term appears on the right hand side of the momentum equations, but does not produce work in the energy equation is a loss mechanism that result in a loss of total pressure and generation of entropy
Diffusion and curvature losses
Viscous losses can be modeled as a contribution of 3 factors:
Conventional skin friction loss
Diffusion loss
Curvature loss
There 3 contributions can be lumped together into the shear stress term τ, which affects both the meridional and tangential momentum components. For some reason, Aungier puts all the loss due to diffusion and curvature into the meridional component. I do not think that there is a strong physical basis for this, and I believe tat this will cause distortions in the exit flow angle.
In order to understand how to model the losses (i.e., loss of total stagnation pressure) we first need to derive the equation for stagnation pressure. Multiplying the momentum equations by their respective velocity components
That is, loss of stagnation pressure is caused by forces that do not contribute to work in the energy equation. τ in this case. Note that the factor 2/b corresponds to the hydrauling diameter and the 1/cosα corresponds to the scale factor between the meridional coordinate m and the length of the streamline s (that is, the losses are proportional to the length of the streamlines, which makes physical sense)
The loss can be modeled as a summ of several contributions and expressed in terms of the pressure loss coefficient
Y=p0−pΔp0
or expressed in differential form
dmdp0=(p0−p)dmdY
where Aungier recommends the following split of the loss
And YC the loss due passage curvature, which si given by:
dmdYF=b226bκ
where κ is the local curvature of the streamline
The term I, as suggested by Aungier describes the additional losses due to diffusion (flow deceleration due to area increase) and curvature of the streamwise channel.
Thoughts for the paper
If there is time, we can explain the model for the 90--degree bend including first the flow equations explaining that they where modified for p-h function calls (in agremeent with look-up tables)
Then explain the geometry.
Simply phi=90 and m=r for the case of radial
for the 90-degree bent we have to define the equations of the NURBS
explain the geometric construction of the shape
3-degree bs-pline for the midline defined by parameters
Also a b-spline (linear) for the thickness distribution.
Explain the equations for the construction and arclength/coordinate reparametrization in an appendix?